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THE MOTION OF A SATELLITE OF THE MOON* 


Harry Lass 

Jet Propulsion Laboratory 
California Institute of Technology 
Pasadena, California 


L INTRODUCTION 


The motion of a satellite of the Moon depends on the potential field due to the Moon as well as the 
gravitational effects of the Earth and Sun. If one chooses a frame of reference attached to the Moon, it can be 
shown that the force field resulting from the Sun can be neglected when compared with the perturbing field of the 
Moon resulting from its oblateness. The effect of the Earth’s field on the satellite is of the same order of 
magnitude as the Moon’s perturbing field and must be included in an analysis of the motion of a satellite of the 
Moon. We will assume that the distance between Earth and Moon remains constant, and we will consider satellite 
orbits of small eccentricity. It will be shown that a nearly circular polar orbit will digress less than 1 deg from 
a polar orbit and that the change in eccentricity is less than a factor of e in one year. 


II. COMPARISON OF THE EFFECTS OF THE EARTH AND SUN 
ON THE MOTION OF A SATELLITE OF THE MOON 


Let m, M, M, Mq be the masses of the satellite, Moon, Earth, and Sun, respectively. We assume that a 
coordinate system centered at the Sun is an inertial frame of reference. From Sketch 1, the motion of the 
satellite is given by 


d 2 r dtp d 2 s 

= m m 

dt 2 dt 2 dt 2 


This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California 
Institute of Technology, under Contract No. NASw-6, sponsored by the National Aeronautics and Space Administration. 
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From 


d 2 p 

m -Jp = F m( M 0^ + F mM + F m(M) 


M — * = f m(M 0 ) + F M (W) 

dr 

with F m (Mq) the force of the Sun on the Moon, etc. (neglecting the force of the satellite on the Moon), Eq. (l) 
becomes 


( 2 ) 


d 2 r _ 

m — = F m^ + 

dt 2 


F m(M 0 ) - ¥ F m(M 0 ) 


F mW ~ ¥ F M<Jf) 


(3) 


CmMn G mM n 


F m (M) + ( - 


P + 


s + 


GmM GmM 

(v + r) + v 


|v+ r| 


The magnitude of the middle term of Eq. (3) can be approximated by GotMq| r |/| s| 3 for |r|« | s|; the 
magnitude of the last term of Eq. (3) can be approximated by GmM\ r\/ \ v| 3 for | r | « | v|. The ratio of these 
terms is 


/Mol 

r 

\»l 

' It 


« (330,000) 

Thus, we are justified in omitting the effect of the Sun. 


240,000 > 
93,000,000 . 


- • 10' 2 

2 


(4) 



Sketch 1 
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III. LAGRANGE’S EQUATIONS OF MOTION 


Let the Lagrangian, L = T - V, of a system of particles be given by L = L(jc x\ x^, x^, ... , t). 

The extremalization of J ( L dt leads to Lagrange’s equations of motion 


d 

{4L\ 

dL 

_ _ Q 

It 

W 

dx i 


1, 2, ••• ,n 


(5) 


We will be interested in transformations of the type 


x* = xHy 1 , 


y", y" +1 » — ,y n+m ) 


* = 1 , 2 , •••,/» 


( 6 ) 


such that the Jacobian \dx l /dy?\, h j = 1, 2, 


n, does not vanish. It is a simple matter to show that 


d 

dt 



dL 

<?y* 



dL dx a 
dx a <?y* 


i = 1, 2, •••, n + m (7) 


with L (y\ •••,y n+m ,y 1 » y n+m , t) = L(x\ ••• , x", x 1 , t). The index, a, of Eq. (7) is summed from 

1 to n. From Eqs. (5), it follows that 


d 

dt 



dL 

— : = 0 

dy* 


i = 1, 2, ••• ,n + m 


( 8 ) 


Moreover, if \dx l /dy?\ ^ 0, i 9 j = 1,2, ••• ,n, one quickly deduces that the vanishing of d/dt(dL/dy *) - dL/dy l 
for i = 1, 2, ••• ,n, yields the set of Eq. (5). Thus y n + \ ••• , y n+m can be chosen in any arbitrary manner, and in 
particular, one can adjoin to Eq. (8) m differential equations involving y 1 , •**,y n+m , provided these new equations 
are not inconsistent with d/dt(dL/dy l ) - dL/dy l = 0, i = 1, 2, ••• , n. We will make use of this result in a sub- 
sequent analysis. 


Jet Propulsion Laboratory 


Technical Release No. 34-56 


IV. THE AVERAGING PROCESS OF KRYLOFF-BOGOLIUBIOFF 


Let x(\p) satisfy the equation 


dx 

Jj, 


£ f {x, sin </», cos t//), 


€« 1 


(9) 


One replaces f(x, sin p, cos </») by its average value, considering * as a parameter, to obtain 


j e f* 

= j sin 1 9, cos 9)dd = eF(x) (10) 

dip 2rr Jq 

We are justified in replacing Eq. (9) with Eq. (10) as x is slowly varying and, hence, remains essen- 
tially constant as granges over the interval (*//, ip + 27r). An integration of Eq. (10) yields x = x(\p, e). This 
method of obtaining an approximate solution of a nonlinear differential equation can be extended to a system of 
differential equations of the type given by Eq. (9). 


4 


Jet Propulsion Laboratory 


Technical Release No. 34-56 


V. QUASI-EULERIAN COORDINATES 

A description of the quasi-Eulerian coordinates associated with a moving point P is given in Sketch 2. 
The x-axis points towards the Earth, the y-axis points in the direction of the orbit of the Moon, and the z-axis is 
the polar axis (axis of rotation of the Moon). The X-Y-Z frame is an inertial frame of reference coinciding with 
the x-y-z frame at t = 0. The unit vector N lies in the x-y plane, and M is a unit vector normal to N in the plane 
formed by N and the position vector r. We define H such that H - N x M. 
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The coordinate transformation between the spherical coordinates and the quasi*Eulerian coordinates 

(a, /?,<//) can be obtained as follows: 


N 

H 


M = HxN 


u 


r 


cos a i + sin a / 

sin /? sin a i - sin /3 cos a / + cos /3 fc 
- sin a cos yS / + cos a cos yS / + sin yS k 

( 11 ) 

cos 0 N + sin 0 M 

(cos 0 cos a - sin 0 sin a cos yS) / + (cos 0 sin a + sin 0 cos a cos yS) / + sin 0 sin /3 k 
sin 8 cos 0 i + sin 8 sin 0 / + cos 0 k 


so that 


sin 8 cos 0 = cos 0 cos a - sin 0 sin a cos yS 

sin 0 sin 0 = cos 0 sin a+ sin 0 cos a cos yS (12) 

cos 0 = sin 0 sin yS 


yielding r = r, 8 = <9 (0, /3), 0 = 0(a,/3,0). This is the type of transformation of coordinates discussed in Sec. 
III. 

The angular velocity of the H, N, M frame of reference is given by 



- cos yS(a + Q ) K + yS N + sin yS ( a + fl ) M 


(13) 


The velocity of the point P is 


dr Dr 

v = = + o) x r 

d£ d£ 


(14) 
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with r - r cos 0 N + r sin 0 M. Thus 


v — r [ sin 0/3 — cos 0 sin y3 ( & + Q )] H + [r cos 0 — r sin 0 0 — r sin 0 cos P (i+n)] N 

(15) 

+ [r sin 0+ r cos 0 0+ r cos 0 cos /3(a+ fi)] M 


which yields the kinetic energy per unit mass given by 


T = — [r 2 + r 2 0 2 + r 2 sin 2 0 /3 2 + r 2 (1 - sin 2 0 sin 2 J3) (a + fl) 2 
- 2 r 2 sin 0 cos 0 sin /3 /3 ( a + ("1 ) + 2 cos /3 0 ( a + 0 )] 


(16) 


After Lagrange's equations of motion are set up, we will impose the constraint that v be perpendicular 
to H, yielding 


sin 0 /3 - cos 0 sin /? ( a + fi ) = 0 


(17) 


Neglecting higher order terms, it can be shown that the Moon's potential is given by (see Ref. 1) 


, CM f* R 3 J ( l . 2 . 2 \ 

1 - (tVu'" 7 


(18) 




[ cos 2a (cos 2 0 - sin 2 0 cos 2 /3) - sin 2 a sin 20 cos y3] ; 


The principal moments of inertia of the Moon about the x-y-z axes, respectively, ore A, B, C. M is the 
mass of the Moon, and R is the radius of the Moon in the direction of the orbit of the Moon (y - axis). The 
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numerical values of / and K are obtained from Ref. 2., although slightly more accurate values can be obtained 
from Ref. 1. The values of / and K are 


/ = 


3 

MR 2 



0.00034 


K = (B - A) m 0.00035 ± 0.000015 

4MB 2 


(19) 


One must also obtain the potential due 
term of Eq . (3) can be written (for a unit mass) 


to the Earth. The x-axis points toward the Earth, thus, the third 
as 


F = 



(x - D) i + y j + z k 
[(* - D ) 2 + y 2 + z 2 ] 2 



( 20 ) 


For * 2 + y 2 + z 2 « D 2 , with D the fixed distance between Earth and Moon, Eq. (20) can be replaced by 


F 


GM 

D 3 


(2 x i-y i 



V 


GM 

2D 3 


( - 3x 2 + r 2 ) 


( 21 ) 


Thus the potential due to the Earth is 

V 2 = — — (r 2 - 3x 2 ) = r 2 (1 - 3 sin 2 6 cos 2 4>) 

2D 3 2D 3 


GM 

2D 3 


r 2 [ 1 - 3(cos 1 // cos a - sin sin a cos /3) 2 ] 


( 22 ) 
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from Eq. (12). The total potential is V = Ej + as given by Eqs. (18) and (22). Comparing R 3 K/r 3 with 
(GM/2D 3 )r 2 for a satellite near the Moon yields 


R 3 K 2D 3 
R 3 GMR 2 


2KD 3 

GMR 2 


2 K D 3 


2 • 3 • 10 " 5 (240, 000) 3 


GM 

R 2 


M 


R 4 


32 

6 - 5280 


= 10 


(23) 


(80) (1100) 4 


Thus, the K-term of V ^ may be slightly more important than the V 2 term in determining the motion of the 


satellite. 
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VI. EQUATIONS OF MOTION 


For / = K = 0, Cl = 0, D = «>, one has Keplerian motion in the M, N plane, provided that H is chosen 
normal to v. For this case, j3 = constant, a= constant, u = 1/r = GM/h^+A sin 0 + B cos 0, with /l, 6, h constants 
of the motion. Thus, for /, K « 1, fi/ 0 « 1, R/D « 1, one assumes that /?, a, A, /l, B will be slowly varying 
quantities. Thus, we will neglect the terms (a+fi) 2 , /3 2 , (a+fi) ft, fi 2 , h a, etc., in Lagrange’s equations of 
motion given by 


d 

dt 



dL 

du 


a = A r 


(24) 


The equation for u = a is a consequence of the Eqs. (24), and will be replaced by Eq. (17). 


The equations of motion are, respectively, 


— [ r 2 sin 2 0/3 - r 2 sin 0 cos 0 sin /3 ( a + fi )] 
dt 


9 „ • , . .. v 2GMR 2 J . 2 f . _ 

- sin p 0 ( a + fi ) sin 0 sin p cos p 


GM/? 2 K o \ • n a . / . /p i 

+ [2 cos 2 a sin^ 0 sin p cos p + sin 2 a. sin 20 sin pj 


3GMr 2 , m f 

+ (cos 0 cos a - sin 0 sin a cos p) sin 0 sin a sin p 

D 3 


(25) 


— [r 2 0 + r 2 cos /3 (a + fi)] 
dt 


2GM/? 2 / . , , . 2 * 2GMR 2 K 

sin 0 cos 0 sin* p 

r 3 r 3 


[ cos 2 a (1 + cos 2 /?) sin 0 cos 0 + sin 2 a cos /3 cos 20] (26) 


3 CMr 2 
D 3 ” 


(cos 0 cos a — sin 0 sin a cos /3) (sin 0 cos a + cos 0 sin a 
cos /3) 
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d 2 r 
it 2 


r 0 2 + 2r cos /30 (a+ fi) - 


GM 


3CMfl 2 / / 1 
\~3 


l 2 <ll ! 


3GMft 2 K 

_ T 4 


[ cos2 a ( cos 2 0 - sin 2 < p cos 2 /3) - sin2 a sin2 0 cos /3] 


(27) 


GMr 


[ 1 - 3 (cos 0 cos a - sin 0 sin a cos /3) 2 ] 


We note that, if /3= 0 at t = 0, then fi= 0 satisfies Eq. (25) and the initial condition, so that an equatorial 
orbit remains equatorial. For (3 4 0, n, making use of Eq. (17) along with r 2 0 = h 4 constant, Eq. (25) becomes 


da 

d\)j 


fir 2 2GMR 2 J . „ , 0 GMR 2 K , , , 

— sin z 0 cos p + (2 cos 2a sin z 0 cos p + sin 2a sin2 0) 

A A2r 


(28) 


3GMr 4 
+ A 2 D 3 


(cos 0 cos a - sin 0 sin a cos /3) sin 0 sin a 


Applying the method of Sec. IV with 


1 GM a . 1 o , CM ( 1 Ah 2 Bh 2 ) 

r h 2 * A 2 \ CM GM I 


Ah 2 


Bh‘ 


CM 


1 sin 0 cos 0 ] for small eccentricity, 

GM GM 


Ah 2 

GM 


« 1 , 


Bh 2 

GM 


« 1. 


one obtains 
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/da\ f H 3 G 2 U 2 R 2 1 G*M 2 R 2 K 3 MA 6 

\ = _ - cos /3 + cos 2a cos y3 — - — - — - cos yS Bin a 

^4>/ C2M2 A 4 *4 2 G 3 M 4 Z) 3 


(29) 


OA 3 C 2 M 2 R 2 


G 2 M 2 A 4 


i yS [ / - K cos 2a + 


SA/A 1 ® 


2G 5 M 6 R 2 D 3 


sin 2 a 1 


From Eqs. (17) and (28), one obtains 


d± 3 

dxjj 


cos yp . _ 

sin p 

sin \p 



2GMR 2 / 

sin yp cos \p cos p sin p 

h 2 r 


2Gm 2 K 

+ 

h 2 r 


( cos 2a sin \p cos >/» sin y0 cos y8 + sin2 a sin (3 cos 2 'Z') 


(30) 


ZGMi 4 

A 2 D 3 


(cos <// cos a — sin <// sin a cos y3) sin a sin [3 cos </» 


Applying the averaging process of Sec. IV yields 


\ 



G 2 M 2 /? 2 / 3MA 10 

K + 

A 4 \ 4G 5 M 6 « 2 D 3 


sin 2 a sin /3 


(31) 
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In order to determine the change in the angular momentum, h = «/» , we return to Eq. (26). Thus, 


dh , a d da nr 2 \ 2GMR 2 ] „ 

= — A cos p I + sin ip cos \p sin z J3 

d\p dip V d\p h ) hr 


2GMR 2 K 

Yr 


[ cos 2a (1 + cos^ /3) sin i p cos i/<+ sin 2a cos /3 cos 2 <J>] 


(32) 


3C/Wr 4 

hD 3 


(cos </z cos a - sin i/z sin a cos /3) (sin «/» cos a + cos t/z sin a cos /3) 


neglecting the terms involving /3 ( a + fi ), A (a + fi ), Making use of Eq. (28) to eliminate d/ dip ( da/dip + fl r 2 IK) 
yields 


dh 4GM/? 2 /u , , , 2GMR 2 J du 2GMR 2 Ku ; 

= cos 2 p sini// cosi// + sin ^ cos * p cosplcos 2 a sin 2ip cosp + sin 2 a cos 2ip) 


d\p h 


h dip 


GMR 2 K du. 


h 

d4> 

12GM 

du 

AD 3 u 5 


2GMR 2 K . 

h 


3GM 

(cos 


3 GM 

AZ)V 


[ sin acosa(cos 2 ip -sin 2 ip) - sin 2 acos/3sin2(//] cos/3 


m f 2GMR 2 J 9 _ 

cos p (cos ip cos a - sin ip sin a cos p) sin ip sin a u sin ip cos ip sin z p (33) 


u [ cos 2 a (1 + cos 2 J3) sin ip cos i// + sin 2a cos /3 cos 2 ip ] 


hD^u^ 


(cos ip cos a - sin ip sin a cos f3 ) (sin ip cos a + cos ip sin a cos /3 ) 
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neglecting the terms involving J ( dh/d\p), K (dk/d\p). Now u = GM/A 2 + A sin \fi + B cos i p, with A, A, B constants 
of the motion for / = K = 0, Q = 0, D = », so that v - du/d\p= A cos i// — S sin \p. We define v - du/d\p ~ 

A cos df - B sin ^// in a subsequent analysis. Applying the averaging process of Sec. IV yields 



(34) 


Thus, for nearly circular orbits, the angular momentum h = r 2 \b remains essentially constant. We are now 
in a position to find a relationship between a and jB from Eqs. (29) and (31). Removing the < > signs one obtains 



3 Mh 10 

4G 5 M 6 /? 2 Z) 3 


sin 2 a sin /3 d a + 


OA 7 

+ cos p 

g 4 m 4 r 2 


/ - K cos 2a + 


m 10 

2G 5 M 6 /? 2 D 3 



d/3 = 0 


(35) 


with h = constant from Eq. (34). Equation (35) becomes exact upon multiplication by the integrating factor sin /3. 
An integration yields 


20. h, 7 

g 4 m 4 r 2 


cos /3 + J-K cos2a + 


3MA 


10 


2G 5 M 6 /? 2 D 3 


sin 2 a ] sin 2 /3 = constant 


(36) 


Accurate observations of a nearly circular satellite of the moon would yield values of a, /?, and A, from 
which one could obtain values of / and K from Eq. (36). 

In order to investigate the stability of a nearly circular polar orbit, let /3= tt/2 - or, cr<< 1, sin /3 = 1, 
cos /3 * a; thus, from Eqs. (29) and (31) 


da GW_ 

dip h 4 

da QA 3 

dip G 2 M 2 



m 


i° \ 


4G 5 M 6 K 2 D 3 j 


sin 2 a 


( 37 ) 
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An integration, with cr= 0 for <Jj = </»q, yields 


a 


c 4 m 4 r 2 


2 flh 7 



m 


10 \ 


4G S M 6 R 2 D 3 ) 


r 

" / (u 3 \ 

< - COS 

2 ro 0 


\ G 2M2 ) 


+ cos 


2 



Clh 3 
G 2 M 2 



(38) 


so that 


< g 4 m 4 r 2 / m 10 \ 

n h 1 l 4G S M 6 R 2 D 3 J 


(39) 


For a satellite near the Moon’s surface, h 2 = G M ft = R 4 2 , < ~ 277/100 rad/ min, (1= 277/ 43000 


rad/ min 


a 


< 


n 




(430) [(3.5) 10' 5 + 6 • 10' 6 ] = (1.5) 10" 2 rad < 1 deg (40) 


thus, an almost circular polar orbit near the Moon’s surface will deviate at most 1 deg from the polar plane. For 
a near circular polar orbit with r = 2 R, 


< 


430 

7 ? 


[(3.5) 10 5 + 2 5 • 6 • 10 6 ] « 8 • 10 3 rad < 0.5 deg 


(41) 


Since the preceding analysis requires an orbit of small eccentricity, it is necessary to study the radial 
motion of the satellite in order to determine the rate of change of the eccentricity, € = h^/CM \j~A^ + # 2 . We 
return now to Eq. (27). From r 2 dj = h 9 it follows that 
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dr h dr ^ du 

dt r 2 d\p dip 


— = - h 2 u 2 — - A U 2 A A 

d £ 2 </i // 2 f/i// rfi/; 


(42) 


Thus, Eq. (27) becomes 


du 

Ip = 


rfi// 


CM 1 eft 

u + v 

A 2 h dp 


2 cos /3 
hu 



3CMR 2 Ju 2 


h 2 



3GNIR 2 Ku 2 

+ 

h 2 


[ cos 2a (cos 2 p - sin 2 ip cos 2 ft) sin2 a sin2 \p cos 0 ] 


(43) 


H [ 1 - 3(cos ip cos a - sin ip sin a cos m 

h 2 D 3 U 3 


The solution of Eq. (43) for / = K - 0, Q. = 0, D = », is 


u 


ai 

77 


+ /4 sin i// + /? cos 0 


(44) 


v = A cos - H sin < p 
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The parameters A and B are now varied, with the result that 


dA 2GM dh 1 dh 2cos/3 / 2 da A 3GMf? 2 / 2 M . o < ? . 

= sin i p v cos ip hu + U cos i p + u sin <p sin p | cos < p 

dip h 3 dip h dip hu \ dip J A 2 W 


+ u 2 [ cos 2a(cos 2 iA- sin 2 i/ 1 cos 2 /3) - sin2a sin2</' cos /3] cos 0 

A 2 


GM 


A 2 D 3 u 3 


[ 1 - 3 (cos ip cos a - sin ip sin a cos ^l 2 ] cos ip 


(45) 


dB 

dip 


2GM dA , 1 dA , 2 cos/3 ( 2 da \ . 3CM/? 2 / 2 /l . 2 , • 2 a\ • , 

cos ip + v sin ip + Au* + U sin ip u ( sin* ip sin* p | sin i// 


A 3 dip A dip hu \ dip 


A 2 \ 3 


3GM/? 2 ft 

A 2 


u 2 [ cos 2 a (cos 2 ip - sin 2 ip cos 2 /3) - sin 2 a sin 2 i p cos /3] sin ip 


GM 

h 2 D*u 3 


[ 1 - 3 (cos ip cos a - sin ip sin a cos /3) 2 ] 


sin ip 


One now substitutes dh/dip, as given by Eq. (33), and [Au 2 ( da/dip ) + 0 ] , as given by Eq. (28), 
into Eqs. (45), with u = CM/A 2 + A sin ip + B cos ip, du/dip = v = A cos ip - B sin ip. Neglecting higher powers 
of Ah 2 /GM, Bh 2 /GM, and applying the averaging process of Sec. IV, one obtains 
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<S>- 


g 2 m 2 r 2 r 


A 4 


^2 - — sin 2 B + K £-3 cos /5 sin 2 a A + — (3 - 5 cos 2 /3) cos 2 a 


M A 6 

15 „ / « 3\ 

+ 

A sin 2 a cos p + B 1 6 cos^ a ) 

d 3 g 3 m 4 

4 \ 2/ 


(46) 1 



g 2 m 2 r 2 r 

~l 



— sin 2 /B'j A + K 


3 cos /3 sin 2 a B + — (3 — 5 cos 2 / 3 ) cos 2 a A 

2 


Mh 6 

D 3 G 3 M 4 


— B sin 2 a cos /3 + (l + cos 2 a. — 5 sin 2 a cos^ fi) I 
_ 4 2 


The terms in < dA/d\p > , <dE/d\fj> involving / and K arose without the assumption of small eccentricity. 
The presence of the Earth’s field on the satellite accounts for the nonsymmetric terms in Eqs. (46); the assump- 
tion of small eccentricity was introduced to derive these terms. R. E. Roberson (Ref. 3) first obtained the critical 
angle given by sin 2 /3= 4/5 (see Eq. 46). 

Removing the < > signs, Eqs. (46) yield 


dA „ dB 3G 2 M 2 /? 2 „ . . n . g d2 n 

A + B K cos /3 sin 2 a ( A 2 + B 2 ) 

djj d^ A 4 


(47) 


Mh 6 

D 3 G 3 M 4 


15 

sin 2 a 

4 


cos /3 C4 2 - B 2 ) + — (cos 2 a- sin 2 a cos 2 /3) AB\ 

2 


For an orbit which remains nearly polar, 77/2, cos /3= 0, and 


1 


F. Yagi, JPL, assisted in checking the results of this equation. 
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dA _ dB 
+ B 

d + d 'P 15 Mh 6 AB < 15 Mh 6 

s: — QQ s " (X _ 

A 2 + B 2 2D 3 G 3 M 4 A2+B 2 4 D 3 C 3 M4 


d\p 2 D3G3M 4 


(48) 


(A 2 + B 2 ) % < / 15 AfA 6 i// \ 

+ = eXp ^oW/ 


with /4 = Aq 9 B = fi 0 for \fj = 0. 

Now the eccentricity of the orbit is given by e = h^/GM \J A^ + ; thus, 


e - e 0 ex P 


15 Mh 6 ' 

4 D 3 G3m 4 


(49) 


with €q the initial eccentricity. The eccentricity will increase at most, by a factor of e when 15/4 Mh 6 /D 2 G 2 M* i// 
= 1, or i p = 4/15 D 3 G 3 M 4 /MA 3 rad. For an orbit near the Moon’s surface, h 2 ~ GMR, and 0 = 4/15 (M/4/) (D/f?) 3 
rads 1/3 • 10 3 rad. The time, T, associated with this increase of eccentricity is 


10 5 100 , 

r = yr ~ 1 year (50) 

677 60 . 24 . 365 


We have shown that a satellite of the Moon will tend to remain in an almost circular polar orbit for a 
considerable length of time. 
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